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Abstract 

■ 

. Following the recent paper by J. F. van Diejen and H. Puschmann we investigate the 

' behavior of zeros of the finite-gap solutions to the Schrodinger equation. As a result, a 

. new system of particles on a punctured Riemann surface is constructed. It is shown to 

be Hamiltonian and integrable. 

p i' 

(D ! 1 Introduction 

It was noticed in [||] that the dynamics of zeros of n-solitonic solutions to the Schrodinger 
equation with the reflectionless potential is governed by the rational Ruijsenaars- Schneider 
system with the harmonic term (||2|). This result appears to be surprising since the afore- 
mentioned dynamics was described long ago, though in the different form. In it was 
shown that the solution to the Schrodinger equation 

{dl - u{x))'il;{x, E) = Eij{x, E) 

with the finite-gap potential u{x) is a well-defined function on the hyperelliptic curve 

2g+l 

y^ = Rg{E)= \{{E-E,). 
1=1 

The projections Q of zeros of this function onto the E'-plane satisfy the Dubrovin equations 
(i): 



dx UjM^s - Cj) 
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Notice that these equations contain the parameters of the curve. The analog of the Dubrovin 
equations holds also for degenerate hyperelliptic curves (the latter are described by the same 
equation where not all i?j's are distinct) and in particular for fully degenerate hyperelliptic 
curves which can be thought of as a Riemann sphere with n couples of pair-wise identified 
points. As it was shown in [P in the latter case the parameters of the curve can be excluded 
from the system. The modified system then is a system of the second-order differential 
equations written solely in terms of zeros of the corresponding function. It coincides with 
the Ruijsenaars- Schneider system and therefore is Hamiltonian, the expressions for the 
parameters of the curve being the integrals of motion. 

In this paper we exploit the algebraic-geometrical approach developed in to apply 
these ideas to the case of the potentials, coming from the hyperelliptic curves with arbi- 
trary degree of degeneracy. The dynamics of zeros of the corresponding solutions to the 
Schrodinger equation is described by the system, which is shown to be Hamiltonian and 
completely integrable, the angle-type variables being the analogs of the components of the 
Abel map. 

In the second section we study the simplest possible case, when a hyperelliptic curve is 
in fact an elliptic curve with several self-intersections. In the third section we generalize the 
results obtained to the case of arbitrary hyperelliptic curve. 

2 Genus one case 

We start with some basic facts from the finite-gap theory. Consider an elliptic curve T, 
given by the equation 

= E^^ - g2E - gs. (1) 

It's compactified at infinity by one point which we denote by oo. The only (up to 

dE 

multiplication by constant) holomorphic differential on F has the following form: uj^ = 

It defines the map from T to the torus T = C/Z[2wi, 2u;2], where 2wi and 2uj2 are o- and 
6-periods of w'*, respectively. This map, given by 




and known as the Abel map, allows us to identify F and F. 

Corresponding to the torus F are the standard Weierstrass functions 

t \ \ r( \ \ ^'(^1^1' ^2) . I s I ^ 

a[z\LJl,i02)., Q{z\iOl,L02) = —r-, r, P{z\^l^^2) = -Q {z\iOl,UJ2) 

a{z\u}i,uj2) 

(see for reference). The function (t{z) has the following properties: 

i) in the neighborhood of zero (t{z) = z + 0{z^); 

ii) a{z + 2ujj) = e2''j(^+'^^ V(z), where rjj = Ci^j)- 

Notice that p{z) is an elliptic function with the only (double) pole at 2; = and C(z) 
has the simple pole at 2; = and satisfies the following monodromy conditions: 

C{z + 2uji) = C{z) + r?i, C{z + 20J2) = Ciz) + m- 
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The map z (E = p{z), y = p'iz)) is inverse to the Abel map. 
Let us fix n — 1 points ki, . . . , on T. 

Proposition 1 For generic divisor = 71 + . . . + 7„ on the curve T there exists a 
unique function 'tp{x,z\D) satisfying the following conditions: 

1. It's meromorphic on the curve T outside the point z = and has poles of at most 
first order at the points 7j, i = 1, . . . , n. 

2. In the neighborhood of z = it has a form 

3. ll){x,Ki) = ll){x,-Ki). 



Remark. In general the function tpi^^ z\D) is defined on the curve F itself, but here for the 
sake of brevity we use the identification between T and F. 

Proof. The uniqueness of such a function follows immediately from the Riemann-Roch 
Theorem. To show the existence we shall consider the following function 

ibix z\D) - e'^(^)^ U7=i'^{z- Zjjx)) 

The set of conditions il^{x, Hi) = ip{x, —Hi) and the constraint Yll=i Zi{x) = x (the latter 
means that ifj is an elliptic function) form the system of n equations on the functions Zi{x). 
For generic data this system is non-degenerate. Then it has the only solution (up to the 
permutations) and therefore defines the function %Ij{x,z\D) uniquely. 

Corollary 1 The above- constructed function iIj{x.,z\D) is a solution to the Schrodinger 
equation 

{dl + u{x))'4}{x, z) = p{z)'4}{x, z), (3) 
where u{x) = —2^27=1 Pi^ii^)) ^i^) ■ 



Proof. Consider a function iPq{x,z) = {d^ -\- u{x) — p{z))ip{x, z). It's straightforward 
to check that the function tjj + ipQ satisfy all defining properties of the function ijj. The 
uniqueness of ip implies that V'o = 0. 

Theorem 1 The zeros of the function ip{x,z\D) satisfy the following dynamics: 

Zj— 2_^ZjZfc — — — — -, « — i,...,n. 

P{zi) - P{Zk) 

Proof. To obtain these equations one has to divide <^ hy ip{x, z) and compare the residues 
of the both sides of the obtained equation at the points Zi{x). 

Remark. Theorem 1 provides us with a wide class of solutions to system (^) coming from the 
algebraic-geometrical data. The simple "dimensional" argument shows that in fact these 
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are all solutions. We could reverse the whole reasoning starting with the solution to (^) and 
showing that the corresponding elliptic function (^) solves the Shrodinger equation. 

^From now on in this section we shall study system (^). Let us introduce the variables 
= In z'^, i = 1, . . . ,n. In the variables Zi, system (^) has the following form: 



Proposition 2 System ^ is Hamiltonian with respect to the Hamiltonian H = ^ e^^ and 
the 2- form 



i=l 



^ = g dz. A d6 - ^ E ^^^^p:^dz. A dz,. (6) 



The proof is a straightforward calculation. 
Note that 

n 



n 



-Pi 

p'{zj) dzj - p'{zi) dzi 



^dzi Ad^i + ^dzi A 



n 



E dzi A d^i + E A diln{p{zj) - pizi))j = ^ dzj A dpi, 

i=l i^j i=l 



where 



Pi=ii + Y^ ln(p(zj) - p{zi)). 

The algebraic-geometrical construction from the previous section provides us with a 
hint on how the first integrals of system (Q) should look like. The constraints ip{x,K,s) = 
^{x,—Ks) imply the equations 



z' 



y 1 = 

f^l P{l^s) - p{zj) 

which can be rewritten in the following form 

n 

J2<Eipi^j)-pi^s)) = o. 

i=l jy^i 

These considerations motivate 

Theorem 2 The coefficients of the polynomial 

n— 1 n 

L(A|z, z') = E H^{z, z')X' = E 4 UiPiz,) - A) (7) 

k=0 i=l j^i 

are the integrals of motion of system (Qj. 
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Remark. Note that the leading coefficient Hn-i{z, z') of L is equal up to the sign to the 
Hamiltonian H(z,z') of system (Q). 

The statement of the theorem is clear since we know that all solutions are algebraic- 
geometrical. However, we would like to present an independent direct proof. It can be 
found in the Appendix I. 

Let us notice that L(p{zj)) = e^K Using this identity we can rewrite the form oj in the 
following way: 

n n 

uj = ^ dzi A dpi = ^ dzi A din L{p{zi)) = 

i=l 1=1 
n 1 /n— 1 \ n n—l g/ \ 



n-l 



dHk A dHs 



n-l 
s=0 



y 

y ^ dE ■ 

^ / L{E)y{E) I 



i=l 



Si 



where the function y{E) is given by (|T|). 

Thus we have proved the following statement. 

Theorem 3 The variables 



^ , E'dE 



i=l 



an 



d Hg defined by ^ are the action-angle type variables for system (QJ. 



We would like however to rewrite the form uj once again in terms of the zeros of the 
polynomial iv(A|z, z') which we shall denote by j = 1, . . . , n — 1. In order to do this we 
introduce the new variables 



1=1 



dE 



{E-Kj)y{E)- 



, n — 1. 



Let us also introduce the variable 



" " dE 



i=l 



y{E) 



Theorem 4 The above-defined form u admits the following representation 

n-l 

UJ = dxA d{ln H) + Y^ dxj A dKj . (8) 
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Remark. We want to emphasize the fact that the variables {x, Xj; J = 1; • • • , — 1} are the 
degenerate curve analogs of the components of the Abel map. So our Hamiltonian structure 
fits in the general scheme proposed in and developed in 

The proof is a straightforward computation (see Appendix II). 

3 General case 

Consider the hyperelliptic curve T of genus g, g > 1, given by the following equation 

2g+l 

T={q = {y, E) y2 = r^^e) = X{{E - Ei)}. 

1=1 

It's compactified at infinity by one point which we denote by oo. The curve F is a 2-sheeted 
branched covering over the complex plane of the variable E. 

Corresponding to the curve F is the following system of differential equations 



Theorem 5 System ^ is Hamiltonian with respect to the 2-f( 



orm 



n 



. A dpj , 



n 

where pj = Cj + ^ ^^{Ck — Cj ) ? o.iT-d the Hamiltonian H = ^ — ^ 



HsiCC); s = 0, . . . , n — 1 of the polynomial 

n—l n A/ 

L(A) = Y: HsiC, C')A^ = E n (C/c - A) 

k = j = l RgiCj) k^j 

are the first integrals of system 1^. Along with the functions 



Coefficients 



n ^j. 



E'dE 

s = 0, . . . , n — 1 



fr{J L{E)y{Ey 
they form the set of action-angle type variables for system Ij^. 

Remark. The level sets of the first integrals are not compact. Thus, there is no canonical 
choice of action-angle type variables for system (P). 

Proof. Let a function Pg{x) be a solution to the differential equation 

dpg 

dx 
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We define new variables Zj, j = 1, ... ,n hy the conditions (j = pg{zj). In terms of these 
variables system (^) has the following form 

The rest of the proof is parallel to the genus one case considered above. 

Remark. Since pg{z) is a well-defined local parameter on the curve F outside infinity, it 
also follows that in fact system (P) describes the motion of n particles on F \ cxo. 

Our next goal is to show that the dynamics of zeros of solutions to the Schrodinger 
equation associated with the curve V is described by the system (^) (the number n of zeros 
being bigger than the genus g of F). 

The way of constructing these solutions is standard in the theory of finite-gap operators. 

Proposition 3 Let's choose a divisor R = ki + . . . + Kn-g on the E -plane. For generic 
divisor D = 71+. . .+7n on the curve F there exists a unique function ^p{x, Q\D, R) satisfying 
the following conditions: 

1. It's meromorphic on the curve F outside the point 00 and has poles of at most first 
order at the points ji, i = 1, . . . ,n. 

2. In the neighborhood of 00 it has a form 

i;{x,Q) = e^^''' (^l + pUx)E-'^'y 

3. 'ilj{x,Kf) = il){x,Hi~), where are preimages of the point Ki under the projection 
onto the E-plane. 

Proof. The uniqueness of such a function (provided it exists) follows immediately from the 
Riemann-Roch Theorem. To show the existence we write this function down explicitly in 
terms of the Riemann 0-function. 

The Riemann 0-function, associated with an algebraic curve F of genus g is an entire 
function of g complex variables z = (zi, . . . , Zg), and is defined by its Fourier expansion 



{zi,...,Zg) = J2 



2iTi(m,z)+TTi{Bm,m) 

m£Z3 



where B = Bij is a matrix of 6-periods, Bij = loj, of normalized holomorphic differen- 
tials ujj{P) on F: Ui = 6ij. Here Oj, bi is a basis of cycles on F with the canonical matrix 
of intersections: Oj o aj = bi obj = 0, aj o bj = 6ij. 

The ^-function has the following monodromy properties with respect to the lattice B, 
spanned by the basis vectors G and the vectors Bj G with coordinates Bij: 

0{z + l) = e{z), e{z + BI) = eicp[-iTT{Bl, I) - 2i7r{l, z)] e{z) 

where / is an integer vector, ^ € Z^. The complex torus J(F) = /B is called the Jacobian 
variety of the algebraic curve F. The vector A{Q) with coordinates 

rQ 

J 10 
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defines the so-called Abel transform: T i-^ </(r)- 

According to the Riemann - Roch theorem, if the divisors D = 71 + . . . + 7„ and TZ = 
+ . . . + , where K'^_g^i = og , are in the general position then there exists a unique 

meromorphic function ra{Q) such that D is its poles' divisor and ra{K^) = S^fs- It can be 

written in the form: 

where 

ff-i 3-1 

s=l s=l 

n n-g+1 

Za = Zo- A{R^), Zo = - /C - ^ ^(7,) + £ ^(iic.), 

s=l a=l 

where /C is the vector of Riemann constants. 

Let dQ be a unique meromorphic differential on T, which is holomorphic outside 00, 
where it has double pole, and is normalized by conditions 

dn = o. 

It defines a vector V with the coordinates 



Define functions ipa{x,Q\D, R) by 

e{A{Q) + xv + z^) e{Zo 



V'a = ra{Q)- 



e{A{Q) + Zc,) e{{xv + z^ 

Let a vector c^, a = 1, . . . ,n — g be a solution to the system 



^ exp ( X I ciri I 
0) V J 



n-g 

X] Calka{Kp) + 1pn+g-l{l^p) = Cp, l3 = 1, . . . ,n - g. 

a=l 

The function 

i;{x, Q\D, R) = ^n-g+i{x, Q\D, i?) + ^ c^M^, Q\D, R) 
satisfies conditions 1-3. 



n-g 
a=l 



Corollary 2 The above-constructed function tp{x,Q\D, R) is a solution to the Schrddinger 
equation 

dl + u{x))^{x,Q)=E^{x,Q), 

where u{x) = —2dx^i{x). 
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The function il^{x, Q) has n zeros on the curve F, whose projections Ci, . . . , satisfy 
the set of Dubrovin equations which in this case read as 



These equations can be rewritten in the matrix form 

Mv = e, (10) 



where 



C' [ (i ^ 3 ^ 9, 



Theorem 6 Projections Cji^), j = o/ zeros o/ the function 'iIj{x,Q\R, D), con- 

structed in the Proposition onto the E-plane satisfy the restriction of system (j^ on the 
joint level set of its g first integrals: 

Hn-k = Skg, k = l,...,g. (11) 
Proof. Consider a function 'ip{x,Q\R, D). Projections Q of its zeros onto E'-plane satisfy 



matrix equation (plj|). The first g equations ensure that restrictions (11) are satisfied. 
Therefore the polynomial L(A|(^, C,') is of degree n — g and the last n — g equations in system 
( [lO|) state that the points ki, . . . , Kn-g are its zeros. The coefficients of L then are time- 
independent, i.e. Hj{x) = 0, j = 1, . . . , n. The last set of equations is equivalent to system 

(ED- 

Remark. In fact using the construction from Proposition ^ we can obtain almost all solutions 
to the restricted system. Namely, for general solution (j, j = 1, . . . ,nto system (P) satisfying 
conditions (p^) we can define algebraic-geometrical data choosing R to be the zero-divisor 



of the polynomial L and D to be (^^i?g(Ci(0)), Ci(0)j + . . . + (^^i?<,(C„(0)), Cn(0) 

Appendix I 

This appendix contains the proof of Theorem ^. 
Let us consider the polynomial 

n n—1 

L(A|z, z') = ^ 4 Uipiz,) - A) = ^ Hk{z, z')XK 

j=l ijLj k=0 

The explicit formulae for the coefficients are 



|J|=n-fc-l jeJ \kfj 
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where summation is taken over all subsets J C {1, . . . , n} of cardinality n — k — 1. 

We are going to show that the functions Hk are time- independent, i.e. dHk/dx = 0. 
Indeed, 



dx 



P{Zs 



k4J 



J jGJ 



J j&J 



\k(^J s&J 



En* 



EE 



(Zs) p{zk) - p{Zs) 



44 



E a{J,k,s), 
J,k(J,seJ 



where 



a{J,k,s) = Yl p{zj) 



p'{Zs) ^ p'jZk) + p'jZs) 
p{Zs) p{zk) - p{Zs) 



I / 

Zj^Zg. 



Let us consider the involution on the set of triples {J, /c ^ J, s G J} which maps {J, /c,s} 
into {J', s, A;}, where J' = J U {A:} \ {s}. 
Now note that 



a( J, A;, s) + a( J', s, /c) 



n p(^i)44[p'(^s) + p'(^fe) + 



p(2;fc) - p{zs) 

and therefore the whole sum J2jk^Jsej'^i'^'^'^) vanishes. 



p{Zs) - p{Zk) 



Appendix II 

Here we present the proof of Theorem ^. 

n-l 

Consider the 2-form uj = E '^Vs A dHs- Recall that Hg = {—ly Han-s-iij^) for s = 

0, . . . , n — 1, where Un-s-iiji) denotes the coefficient of A'' in the polynomial Hi^T/l'^ + ^i)- 
By cr^_g_2('?) we denote the coefficient of A'' in the polynomial Ili^jl-^ + ^i)- Then 

n—l n—2 

u = Y.d^s^dHs = J2'^'P^^i-'^y'^(H'^n-s-iifi)) + i-ir~^dipn-iAdH = 

n-l n-ln-2 U^j 

= E d^s^dH+J2Yl {-'^yHa{_^_^{K) d^s A 

s=0 j=l s=0 



10 



Now let us notice that 



E'dE 



n—2 n—2 n 

Y^{-irHai_,_,{2)d^s = Y.lli-^yH<.,^2r^)d 

.=0 8=0 1=1 ^ HJ^)y[E) 

/n-2 \ 



1=1 



n-2 
s=0 



L{E)y{E) 



n 

dE-J2 / d 

1=1 



s=0 



L{E)y{E) 



1 1 /I 

^'1 



dE 



V 



dE 



+ E 



dE 



In the same way one can show that 



dKj. 



n—l n * Y 

'^{-iyan-s-i{K)d(ps = '^d / 

s=0 1=1 



dE 



p{^i) 



+ E 



Hy{E) J my{E) 



dE 



dH. 



Plugging these two formulae into (12) we obtain 

= EdiT. l^^ \^_7A..f^., ]^dK,+d(±dJ 



n-l 



AdH 



= ^ dxj A dKj + dxA d(ln H) . 
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